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Abstract — The operations managers of Mega Petroleum 
Stations were not able to determine the best number of servers 
that can serve arriving customers at various demand periods 
which affects their queue performance. This study was 
conducted at NNPC Mega Petroleum Station in Enugu, Nigeria 
with the aim of addressing the identified problem. Experimental 
observations were conducted simultaneously on the referenced 
service facilities in order to collect the daily arrival rates of 
customers. Arrival rates and combined service rates were 
however, collected at every 15 minutes interval and at peak 
demand periods of the day. From the results of the queuing 
evaluations, it was discovered that at an average of 6 servers 
being used, with combined service rate () of 1.5320 

cars/minutes and average customer arrival rate (X) of 1.5153 
cars/minutes, for NNPC mega petroleum station Enugu gave a 
system utilization (P) of 0.9892 which gave a percentage system 
utilization of 98.9%. It was discovered that the service systems 
were being over utilized at almost 100% which resulted to the 
longer waiting time of customers at the service facilities. The 
result from the Queue Evaluation Environment showed that 8 
servers gave the best system utilization value of 74.2% which 
reduced the customers waiting times (Ws) by 92.7%. The 
expected probability of system idleness for the case study is 
negligible at 8 server utilization. The Queue Evaluation 
Environment was later adopted in developing a Decision 
Support System for the referenced service facilities. The 
Decision Support System was finally recommended to guide the 
operations manager in determining the best number of servers 
to engage at various demand periods for PMS refill only. 

Index Terms — Queuing, NNPC Mega Station, Server, PMS, 
Decision support system, waiting time, system utilization, 
Arrival rate and Service rate. 


I. INTRODUCTION 

Queuing theory is the mathematical study of waiting lines 
[1]. The theory permits the derivation and calculation of 
several performance measures which includes the average 
waiting time in the queue or the system, the expected number 
waiting or receiving service, the probability of encountering 
the system empty, having an available server or having to wait 
a certain time to be served and most importantly the system 
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utilization [2]. As a result of its applications in industries, 
technology, telecommunications networks, information 
technology and management sciences, it has been an 
interesting research area for many researchers active in the 
field. 

The theory of queues was initiated by the Danish 
mathematician Erlang, who in 1909 published “The theory of 
Probabilities and Telephone Conversation”. He observed that 
a telephone system was generally characterized by either (1) 
Poisson input (the number of calls), exponential holding 
(service) time, and multiple channels (servers), or (2) Poisson 
input, constant holding time and a single channel. Erlang was 
also responsible in his later works for the notion of stationary 
equilibrium and for the first consideration of the optimization 
of a queuing system. 

In 1927, Molina published “Application of the Theory of 
Probability to Telephone Trunking Problems”, and one year 
later Thornton Fry printed “Probability and its Engineering 
Uses” which expanded much of Erlang’s earlier work. 
Kendall was the pioneer who viewed and developed queuing 
theory from the perspective of stochastic processes [3]. 
Kleinrock also did some extensive work on the theory of 
queuing systems and their computer applications [4]. The 
work in queuing theory picked up momentum rather slowly in 
its early days, but in 1960’s started to accelerate and there 
have been a great deal of work in the area and its applications 
since then [5]. 


II. Review of Related Work 

In recent times, queuing theory and the diverse areas of its 
applications has grown tremendously. Takagi considered 
queuing phenomena with regard to its applications and 
performance evaluation in computer and communication 
systems [6]. Obamiro Applied Queuing Model in 
Determining the Optimum number of Service Facility needed 
in Nigerian Hospitals. He however achieved this by 
determining some queuing parameters which enabled him to 
improve the performance of the system [7]. Azmat also 
applied queuing theory to determine the sales checkout 
operation in ICA supermarket using a multiple queue multiple 
server model. This was used to obtain efficiency of the models 
in terms of utilization and waiting length, hence increasing the 
number of queues so customers will not have to wait longer 
when servers are too busy. The model contains five (5) 
servers which are checkout sales counters and it helps to 
reduce queue [8]. 
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Yankovic N. and Green L. developed a queuing model to help 
identify nurse staffing levels in hospital clinical units based on 
providing timely responses to patient needs. The model 
represents the crucial interaction between the nurse and bed 
systems and therefore includes the nursing workload due to 
admissions, discharges and transfers, as well as the observed 
impact of nursing availability on bed occupancy levels [9]. 
Mgbemena was able to model the queuing system of some 
banks in Nigeria using regression analysis. In her work, she 
created queuing management software in MATLAB that 
shows at a glance, the behavior of the queuing system and the 
unit that needs attention at any time. The essence was to 
improve the customer service system in Nigerian banks [10]. 
Vasumathi and Dhanavanthan applied Simulation Technique 
in Queuing Model for ATM Facility. The main purpose of 
their study was to develop an efficient procedure for ATM 
queuing problem, which can be daily used by banks to reduce 
the waiting time of customers in the system. In their work, 
they formulated a suitable simulation technique which will 
reduce idle time of servers and waiting time of customers for 
any bank having ATM facility [11]. 

Ahmed S.A and Huda K.T focused on banks lines system, the 
different queuing algorithms that are used in banks to serve 
customers, and the average waiting time. The aim of their 
paper was to build an automatic queuing system for 
organizing and analyzing queue status and take decision of 
which customer to serve in banks. The new queuing 
architecture model can switch between different scheduling 
algorithms according to the testing results and the factor of the 
average waiting time. The main innovation of their work 
concerns the modeling of the average waiting time taken into 
processing, in addition with the process of switching to the 
scheduling algorithm that gives the best average waiting time 
[ 12 ]. 

Chinwuko and Nwosu adopted the single line multi-server 
queuing existing model to analyze the queuing system of First 
Bank Nigeria PLC. In their work, they suggested the need to 
increase the number of servers in order to serve customers 
better in the case study organization [13]. Tabari et al used 
queuing theory to reorganize the optimal number of required 
human resources in an educational institution carried out in 
Iran. Multi-queuing analysis was used to estimate the average 
waiting time, queue lengths, number of servers and service 
rates. The analysis was performed for different numbers of 
staff members. Finally, the result shows that the staff members 
in this department should be reduced [14]. 

Ohaneme et al, proposed the single line multi-server queuing 
system which they simulated using c-programming to be 
adopted at NNPC Mega petroleum station in Awka, Anambra 
State in order to avoid congestion and delay of customers [1]. 
Akpan N.P et al, studied queuing theory and its application in 
waste management authority in LAWMA Igando dump site, 
Lagos state by adopting the M/M/S queuing model. He used 
the Queuing performance measures to estimate the 
inter-arrival, service and waiting time of the queue. His study 
showed that both the service and the inter-arrival time made a 
good fit to Exponential distribution. 

However, this work goes further in evaluating the 
performance of the queuing system, creating a Queue 
Evaluation Environment that gives expected queue 
performance and developing a Decision Support System that 
recommends the best number of servers to use at various 
demand periods [15]. 


The objective of the study is to address the queuing problem 
at NNPC Mega Petroleum Stations by developing a Decision 
Support System that recommends the best number of servers 
needed to be engaged at various demand periods. 

The structure of the studied system is shown in figure 1. The 
structures can be approximated as a single-line multi-server 
queuing systems. . At the NNPC mega petroleum station, 
there are six dispensers i.e. fuel metering pumps (5 1 to 5 & ) 
in the system. Each of the fuel dispensers has two nozzles. 
This means that at full capacity of operation the service 
facility should be considered as a twelve-server system. 


III. Methodology 

The research method used in this work was the quantitative 
research approach. . The single line multi-server queuing 
model was adopted for developing the results of the queue 
performance. This model was adopted because it showed a 
good representation of the model structure of both case 
studies of queuing systems. 

3.1 Method of Data Analysis 

The data generated was first organized and descriptive 
statistics was used to compute the total average arrival rates 
and total average combined service rates for the year. The 
service rates per server of both facilities were established and 
the single line multi server queuing model was coded in 
Microsoft Excel using 2-12 servers (i.e. when M = 2 - 12 
servers) in creating the Queue Evaluation Environment that 
generates the expected queue performance results at the 
respective average arrival rates of customers in the referenced 
service facilities. The flow chat developed of the Queue 
Evaluation Environment is presented in figure 2. The Queue 
Evaluation Environment was later adopted in developing the 
decision support system using the application of Microsoft 
Excel. 


3.2. Models Applied for the Queuing Analysis 

Based on the assumptions of the single line multi-server 
queuing model, the expressions for the performance measures 
which are derived from the analysis of the birth-and-death 
models, [16, 17andl8] are; 

i. The average utilization of the system: 


When m = 6 is 


p= ^ 

When m = 2- 12 is 

p = 

M (p) 


( 1 ) 


( 2 ) 


ii. The probability that there are no customers in the 
system is 

pf pf 

?M- lW i W 


p q = 


ywr-i 

*-‘n=a 


(3) 


+ 


< 4 -£) 


iii. The average number of customers waiting for service. 

7 af 

l = _L±iJ_ p (4) 
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iv. The average number of customers in the system. 

L a = L q + (X/i 0 (5) 

v. The average time a customer spends in line waiting for 

service 

w q = ( 6 ) 

vi. The average time a customer spends in the system. 

w s = \ ( 7 ) 

vii. The average waiting time of a customer on arrival 

not immediately served. 


viii. 


K = 


1 

MpL-2 


( 8 ) 


Probability that an arriving customer must wait 


P 

w w„ 


(9) 


It is seen that these performance measures depend on two 


basic queue parameters, namely; ^and/T Given A and/t, the 
values computed for these measures gives an indication of 
how well the referenced service facilities handle the volume 
of arriving customers. 



Figure 1: Structure of the PMS Dispensary pump system of the studied NNPC mega petroleum station Enugu 
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Figure 2: Development Flow Chart of the Queue Evaluation Environment and Best Server Utilization. 
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IV. Data Analysis (Case Study: NNPC Mega 
Petroleum Station Enugu) 



Figure 3: Total Daily Average Arrival Rate /15 Minutes 


In figure 3, the bar chart shows the total daily average arrival 
rate of customers Per 15 Minutes for the year from Monday - 
Sunday at NNPC mega petroleum station Enugu. From the 
chart, it is observed that Saturday was with the highest arrivals 
which shows that the mega station is being patronage more by 
customers on Saturdays than the rest of the days being the fact 
that Saturday is a work free day for civil servants and most 
public servants so customers buy large quantity of PMS to last 
them for the week day activities as well as weekend travels. 
While Sunday was with the lowest arrivals which shows less 
patronage of customers being the fact that Sunday is a worship 
day for Christians and the Mega Station don’t always open for 
service on that day and also customers most especially civil 
and public servants must have bought large quantity of PMS 
on Saturday and Friday to last them for the week. The week 
days (i.e. Mondays - Fridays) were mostly patronage more by 
commercial transporters 
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1.5153 
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6 


Total average arrival rate for the year 

q-Vt-gn lj. j _ 5^53 cars/minutes 

24 


Total average number of servers being used for the 


Table 1: Weekly Mean Arrival Rate of Customers, 
Weekly Mean Combined Service Rate of Customers and 
Mean Number of Servers Engaged at NNPC Mega 




Weekly 

Mean 

Arrival 

Weekly 
Mean 
Combined 
Service 
Rate Per 
Mins 

Mean 
Number of 
Servers 



Rate Per 
Mins 

Being 

Used (M) 


1st Week 

1.5488 

1.5571 

6 

DEC 

Last 

Week 

0.3893 

0.4226 

4 
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0.9036 

0.9554 

5 
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1.5804 

6 


1st Week 

1.5601 

1.5738 

6 

MAR 

Last 

Week 

1.6857 

1.7077 

6 


E* rw-rafrbr 'i 

year =-'- 

J 24 


= 6 


Total average combined service rate for the year = 

(W-ggmjj _ j 5320 cars/minutes. 

24 

(Note) It is assumed that each server contributes an average 
service rate of —■ ~ ■ cars/minutes. Where M = 6, and 

j't? 

{■ Ye ar) = 1-^20 cars/minutes. This implies that each server 

contributes an average service rate of 0.2553 cars/minutes in 
the service facility 


V. Presentation of Queue Performance Evaluation 
Results 

The results of the performance measures of the queuing 
system are presented below. 

Table 2: Results of the performance evaluation of the 
queuing system with parameters 

A = 1.5153 cars/minutes and fi c = 

1.5320cars/minuteswhen (Af= 6) 
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PvsM 
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y = 5.9354X 1 
R 2 = 1 


Power (P) 


M 


From figure 4, the charts of the queue output results were 
developed using the application of Microsoft Excel and trend 
line was used to test for the best goodness fit in developing the 
relationship that exists best between the queue output results. 
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Figure 4: Queue Evaluation Environment displaying the 


results of the queue performance when A is fixed, 
0.2553 cars/minutes (per server) and M- 2-12 servers 
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Figure 5: Scatter Plot of System Utilization vs. Number of 
Servers 


In figure 5, the scatter plot shows that the best fit between the 
two variables i.e. (P and M) from 2 to 12 servers is a trend line 
power. _ 

P 0 vsM 

y = -3E-05X 6 + 0.0015X 5 -0.0282x 4 + 0.2761X 3 - 



Figure 5: Scatter Plot of Probability System is Empty vs. 
Number of Servers 
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In figure 5,the scatter plot shows that the best fit between the L q = -0.012m 6 + 0.471m 5 - 6.701x 4 + 44.44m 3 - 137.6m 2 + 
two variables i.e. (P 0 and M) from 2 to 12 servers is a 177.5m-71.61 (10) 

nonlinear polynomial function in the sixth order. R 2 = 0.360 


P w vsM 



Figure 6: Scatter Plot of Probability Arrival Must Wait 
vs. Number of Servers 


In figure 6,the scatter plot shows that the best fit between the 
two variables i.e. (P w and M) from 2 to 12 servers is a 
nonlinear polynomial function in the third order. 
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Figure 7:Scatter Plot of Average Number in Line, 
Average Number in System, Average Time in Line, 
Average Time in System and Average Waiting Time vs. 
Number of Servers. 


In figure 7, the scatter plot shows the relationship that exists 
between the queue output variables i.e. (L q , L s , W q , W s and 
W a plotted against M) from 2 to 12 servers. The chart can 
hardly be interpreted and the best fit between the variables is a 
nonlinear polynomial function in the sixth order which 
produced poor R 2 values at 0.360 for L q , L s , W q , W s and 0.356 
for W a . (See equation 10 - 15). The R 2 values are poor and 
thus, the relationship that exists between each of the 
dependent variables i.e. (L q , L s , W q , W s and W a ) and the 
independent variable (M) were poor. This resulted from the 
negative values of L q , L s , W q , W s , and W a from 2 to 5 servers 
(see figure 1). The resultant of these negative values shows 
that system utilization is greater than 1. The scatter plots of 
the variables were later plotted separately from the positive 
result outputs (i.e. from 6 to 12 servers) and the results 
produced excellent R 2 values. (See figure 8-12). 


L s = -0.012m 6 + 0.471m 5 - 6.701m 4 + 44.44m 3 - 137.6m 2 + 
177.5m-65.67 (12) 

R 2 = 0.360 

W q = -0.008m 6 + 0.31 lm 5 - 4.422m 4 + 29.33m 3 - 90.84m 2 + 
117.1m-47.26 (13) 

R 2 = 0.360 

W s = -0.008m 6 + 0.31 lm 5 - 4.422m 4 + 29.33m 3 - 90.84m 2 + 
117.1m-43.34 (14) 

R 2 = 0.360 

W a = -0.008m 6 + 0.31 lm 5 - 4.422m 4 + 29.33m 3 - 90.82m 2 + 
116.6m-42.81 (15) 

R 2 = 0.356 


L q vsM 



Figure 8: Scatter Plot of Average Number in Line vs. 
Number of Servers 

In figure 8, the scatter plot shows that the best fit between the 
two variables i.e. (L q and M) from 6 to 12 servers is a 
nonlinear polynomial function in the sixth order. 


L s vsM 



M 


Figure 9: Scatter Plot of Average Number in System vs. 
Number of Servers 
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In figure 9, the scatter plot shows that the best fit between the 
two variables i.e. (L s and M) from 6 to 12 servers is a 
nonlinear polynomial function in the sixth order. _ 


W q vsM 



Figure 10: Scatter Plot of Average Time in Line vs. 
Number of Servers 

In figure 10, the scatter plot shows that the best fit between the 
two variables i.e. (W q and M) from 6 to 12 servers is a 
nonlinear polynomial function in the sixth order. 


W s vs M 



Figure 11: Scatter Plot of Average Time in System vs. 
Number of Servers 



Figure 13: Scatter Plot of Probability System is Empty vs. 
System Utilization 

In figure 13, the scatter plot shows that the best fit between the 
two variables i.e. (P 0 and P) from 2 server utilization to 12 
server utilization is a nonlinear polynomial function in the 
third order. 


P„vsP 



p 


Figure 14: Scatter Plot of Probability Arrival Must Wait 
vs. System Utilization 


In figure 11, the scatter plot shows that the best fit between the 
two variables i.e. (W s and M) from 6 to 12 servers is a 
nonlinear polynomial function in the sixth order. 



Figure 12: Scatter Plot of Average Waiting Time vs. 
Number of Servers 

In figure 12, the scatter plot shows that the best fit between the 
two variables i.e. (W a and M) from 6 to 12 servers is a 
nonlinear polynomial function in the sixth order. 


In figure 14, the scatter plot shows that the best fit between the 
two variables i.e. (P w and P) from 2 server utilization to 12 
server utilization is a nonlinear polynomial function in the 
fifth order. 


Figure 15: Scatter Plot of Average Number in Line, 
Average Number in System, Average Time in Line, 
Average Time in System and Average Waiting Time vs. 
System Utilization 
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In figure 15, the scatter plot shows the relationship that exists 
between the queue output variables i.e. (L q , L s , W q , W s and 
W a plotted against P) from 2 Server Utilization to 12 Server 
Utilization. The chart can hardly be interpreted and the best fit 
between the variables is a nonlinear polynomial function in 
the sixth order which produced poor R 2 values at 0.430 for L^, 
L s , W q , W s and 0.426 for W a (See equations 13 - 17). The R 2 
values are poor and thus, the relationship that exists between 
each of the dependent variables i.e. (L q , L s , W q , W s and W a ) 
and the independent variable (P) were poor. This resulted 
from the negative values of L q , L s , W q , W s , and W a from 2 to 5 
Server Utilization (see figure 3). The resultant of these 
negative values shows that system utilization is greater than 1. 
The scatter plots of the variables were later plotted separately 
from the positive result outputs (i.e. from 6 servers’ utilization 
to 12 servers’ utilization) and the result produced excellent R 2 
values. (See figure 16 - 20) 

L q =795.7p 6 - 7176.p 5 + 25352p 4 - 44779p 3 + 41488p 2 - 
19002p + 3360. (16) 

R 2 = 0.430 

L s = 795.7p 6 - 7176.p 5 + 25352p 4 - 44779p 3 + 41488p 2 - 
19002p + 3366. (17) 

R 2 = 0.430 

W q =525.Ip 6 - 4736.p 5 + 1673lp 4 - 2955lp 3 + 27379p 2 - 
12540p + 2217. (18) 

R 2 = 0.430 

W s =525.lp 6 - 4736.p 5 + 1673lp 4 - 2955lp 3 + 27379p 2 - 
12540p + 2221. (19) 

R 2 = 0.430 

W a = 525.4p 6 - 4739.p 5 + 16743p 4 - 29574p 3 + 27401p 2 - 
12547p + 2218. (20) 

R 2 = 0.426 


C VSP 



Figure 16: Scatter Plot of Average Number in Line vs. 
System Utilization 

In figure 16, the scatter plot shows that the best fit between the 
two variables i.e. (L q and P) from 6 server utilization to 12 
server utilization is a nonlinear polynomial function in the 
sixth order. 



Figure 17: Scatter Plot of Average Number in System vs. 
System Utilization 


In figure 17, the scatter plot shows that the best fit between the 
two variables i.e. (L s and P) from 6 server utilization to 12 
server utilization is a nonlinear polynomial function in the 
sixth order. 


W q VS P 



Figure 18: Scatter Plot of Average Time in Line vs. 
System Utilization 

In figure 18, the scatter plot shows that the best fit between the 
two variables i.e. (W q and P) from 6 server utilization to 12 
server utilization is a nonlinear polynomial function in the 
sixth order. 



Figure 19: Scatter Plot of Average Time in System vs. 
System Utilization 


In figure 19, the scatter plot shows that the best fit between the 
two variables i.e. (W s and P) from 6 server utilization to 12 
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server utilization is a nonlinear polynomial function in the 
sixth order. 


w a vsP 



Figure 20: Scatter Plot of Average Waiting Time vs. 
System Utilization 


In figure 20, the scatter plot shows that the best fit between the 
two variables i.e. (W a and P) from 6 server utilization to 12 
server utilization is a nonlinear polynomial function in the 
sixth order. 


L q vsW q 



y= 1.5151X-0.0004 
R 2 = 1 


-Lq 

-Linear (Lq) 


Figure 21: Scatter Plot of Average Number in Line vs. 
Average Time in Line 


In figure 22, the scatter plot shows that the best fit between the 
two variables i.e. (L s and W s ) from 2 to 12 servers is a linear 
function. 


VI. Development of the Decision Support System for 
Case Study A and B 


From the Queue Evaluation Environment created, the service 
rates (per server) of each of the referenced facilities were 
fixed and arrival rates were simulated using 2-12 servers to 
see the expected queue performance and to determine the best 
number of servers that gives the best system utilization value 
at various arrival rates of customers. The summary result 
outputs were plotted on a chart using the application of 
Microsoft Excel and trend line was used to test for the best 
goodness fit between the dependent variable i.e. Number of 
Servers (M) and the independent variable i.e. Average Arrival 
Rates/Minutes (A). See summary result output and charts 
below. 

Table 3: Summary result output of simulated arrival 
rates of customers/minutes (NNPC Mega Petroleum 
Station Enugu) 

Af Mins Best No. Servers 
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In figure 21, the scatter plot shows that the best fit between the 
two variables i.e. (L q and W q ) from 2 to 12 servers is a linear 
function. 


Ls vs Ws 


y = 1.5153X+7E-05 



Figure 22: Scatter Plot of Average Number in System vs. 
Average Time in System 


M vs A/Mins 


y = 1.1908X 4 - 5.2393X 3 +8.0734x 2 - 0.0822x+1.3308 
R 2 = 0.9902 



Average Arrival Rate/Minutes 


Figure 23: Scatter Plot of Number of Servers (M) vs. 
Average Arrival Rate/Minutes 


From figure 23, the scatter plot shows the number of servers 
plotted against average arrival rate of customers/minutes in 
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NNPC Mega Petroleum Station Enugu. From the chart, it is 
observed that number of servers is expected to be stepping up 
as average arrival rate increases and the best fit between the 
two variables i.e. Number of Servers (M) and Average Arrival 
Rate (X) is a nonlinear polynomial function in the fourth order 
as depicted in the chart. 

VII. Discussion of Results 

From the analysis, table 2 showed the results of the 
performance measures of the queuing system as seen at the 
NNPC mega petroleum station Enugu. From the results, it was 
discovered at an average number of 6 servers being used, with 
average combined service rate(/I c ) of 1.5320 cars/minutes 
and average customer arrival rate (A) of 1.5153 cars/minutes, 
gave a system utilization (P) of 0.9892 which gives a 
percentage system utilization of 98.92%, while the 
probability of the system being empty and the probability of 
waiting gave 0.0002 and 0.9703 respectively, this means that 
when service commences, the system is never idle and a 
customer must wait before receiving service with a 97.03% 
probability. However, the average number of customers in 
line and the average number of customers in system including 
any being served gave 89.1 and 95.045 respectively. 
Furthermore, the average waiting time of customers in line, 
the average waiting time of customers in the system including 
service and the average waiting time of customers on arrival 
not immediately served gave 58.807, 62.724 and 60.606 
minutes respectively. 

The results of table 2 studies showed that the systems were 
heavily utilized at an average of 6 servers because system 
utilization was almost 100%. This resulted to the longer 
waiting time of customers experienced at the service facilities. 
However in respect of this, the service rate per server was 
determined for the studies and a Queue Evaluation 
Environment was created using 2-12 servers to see the 
expected queue performance and to determine the best 
number of servers that gives a good tradeoff between system 
utilization and waiting time at the collected average arrival 
rates of customers in the referenced service facilities. 

The results from the Queue Evaluation Environment showed 
that 8 servers gave the best system utilization values of 0.7419 
which is expected to reduce the respective customers waiting 
times (Ws) by 92.72% for the study establishments. This is 
based on the statement of Egolum, which says that system 
utilization should be greater than 0 but less than 0.8 [19]. 
From the charts of system utilization versus waiting time 
plotted for both case studies, it is observed that there’s no 
significant decrease in waiting time anymore from system 
utilization value of 0.8, which shows that waiting time has 
reached its optimum at the respective best server utilization 
values of 0.7419 of the referenced service facilities. See 
figures (18, 19 and 20). This shows that there will be no need 
of making use of more than 8 servers at the respective average 
arrival rates of customers in the referenced service facilities. 
Also the expected probabilities of system idleness for the 
study is negligible at 8 server utilization because at that point, 
probabilities of system idleness has also reached its optimum 
and it no longer has any effect on the service systems. (See 
figure 13). 

From the study, it was revealed that the probability of system 
being empty increased to optimum as system utilization 


drops; the probability of an arrival waiting drops as system 
utilization reduces; the average number in line and average 
number in system drops to optimum as system utilization 
reduced; the average time in line, average time in system and 
average waiting drops to optimum as system utilization 
reduced; the average number in line increases as average time 
in line increases; the average number in system increases as 
average time in system increases (See figures (13-14), (lb- 
20) and (21 - 22). 

It was also revealed that system utilization drops as number of 
server’s increases; the probability of system being empty 
increased to optimum as number of server’s increases; the 
probability of an arrival waiting reduces as number of server’s 
increases; the average number in line and average number in 
system drops to optimum as number of server’s increased; the 
average time in line, average time in system and average 
waiting time drops to optimum as number of server’s 
increased (See figure 8-12). 

Finally, from figure 23, the models for the decision support 
system were developed using trend line analysis. For NNPC 
Mega Petroleum Station Enugu, the recommended best 
number of server (M) for PMS refill is given by: 

M = 1.190A 4 - 5.239V + 8.073V - 0.082A + 1.330 

( 21 ) 

VIII. Conclusion 

The evaluation of queuing system in an establishment is very 
essential for the betterment of the establishment. Most 
establishments are not aware on the significance of evaluating 
their queue performance. The implication of this is that 
operations managers are not able to determine the best 
number of servers to engage for service at various demand 
periods which affects their queue performance. As it concerns 
the case study establishments, the analysis and evaluation of 
their queuing system showed that both service systems were 
being over utilized which resulted to customers spending 
longer time than necessary before receiving service. 
However, the need of creating a Queue Evaluation 
Environment to find out the number of servers that gives the 
best server utilization at the collected average arrival rates 
became very essential. From the Queue Evaluation 
Environment, using 8 servers at the collected average arrival 
rates of customers in the referenced service facilities gave a 
good tradeoff between system utilization and waiting time 
which is expected to reduce the waiting time of customers’ 
inthe system while server idleness is neglected. In conclusion, 
the Queue Evaluation Environment created and the decision 
support system developed for the case study establishments 
will go every long way in addressing their queuing problems. 
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